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In this paper, we study the three-body decays B^/B^ —)• pcfoi^) 77c7r+7r“ by em¬ 
ploying the perturbative QCD (PQCD) factorization approach. We evaluate the 5-wave res¬ 
onance contributions by using the two-pion distribution amplitude The Breit-Wigner 
formula for the /o(500), /o(1500), and /o(1790) resonances and the Flatte model for the 
/o(980) resonance are adopted to parameterize the time-like scalar form factors Fsico^). We 
also use the D. V. Bugg model to parameterize the /o(500) and compare the relevant theoret¬ 
ical predictions from different models. We found the following results: (a) the PQCD pre¬ 
dictions for the branching ratios are B{B^ —)• ryc/o(500)[7r'''7r“]) = (l.SS^Qg®) x 10“® for 
Breit-Wigner model and B{B^ —)• ryc/o(500)[7r'''7r“]) = (2.31^0 Ig) x 10“® for D. V. Bugg 
model; (b) B{Bs —)> ?7c/o(X)[7r'''7r“]) = (5.02l} Qg) x 10“^ when the contributions from 
fo{X) = (/o(980), /o(1500), /o(1790)) are all taken into account; and (c) The considered 
decays could be measured at the ongoing LHCb experiment, consequently, the formalism of 
two-hadron distribution amplitudes could also be tested by such experiments. 

PACS numbers: 13.20.He, 13.25.Hw, 13.30.Eg 


I. INTRODUCTION 

Several years ago, some three-body hadronie B —)■ 3h^'^ ( h,h' = tt, K) deeays have been 
measured by BaBar and Belle Collaborations [ 1 ] and studied by using the Dalitz-plot-analysis. The 
LHCb Collaboration reported, very recently, their experimental measurements for the branching 
ratios and sizable direct CP asymmetries for some three-body charmless hadronic decays —)■ 
K~, and tt+tt+tt" [2-4], the three-body charmed hadronic decays 

B^ —)■ J/tpTt^Tt~ [5-8] and —)■ [9], or the decays B^ —)■ J/ficjxl), 

D^K~7r~^ [6, 7, 10]. The large localized CP asymmetry observed by LHCb brings new challenges 
to experimentalists and their traditional models to fit data, and also has invoked more theoretical 
studies on how to understand these very interesting three-body B/Bg meson decays. 

On the theory side, the three-body hadronic decays of the heavy B/Bg meson are much more 
complicated to be described theoretically than those two body B/Bg —)■ hih 2 (here hi refer to 
light mesons vr, K, p, etc ) decays. During the past two decades, such two-body hadronic B/Bg 
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meson decays have been studied systematically and successfully by employing various kinds of 
factorizations approaches. The three major factorization approaches are the QCD-improved fac¬ 
torization (QCDF) [11-13], the perturbative QCD (PQCD) factorization approach [14-18] and 
the soft-collinear-effective theory (SCET) [19, 20]. For most B/Bg —)■ hih 2 decay channels, the 
theoretical predictions obtained by using these different factorization approaches agree well with 
each other and also be well consistent with the data within errors. 

For B/Bg three-body hadronic decays, however, they do receive both the resonant and non¬ 
resonant contributions, as well as the possible final state interactions (FSIs), while the relative 
strength of these contributions are varying significantly for different decay modes. They are known 
experimentally to be dominated by the low energy resonances on tttt, KK and Kir channels on 
Dalitz plot, usually analysed by employing the isobar model in which the decay amplitudes are 
parameterized by sums of the Breit-Wigner terms and a background, but without the inclusion 
of the possible contributions from the coupled channels and the three-body effects such as the 
FSIs. In fact, the three-body hadronic B/Bg meson decays have been studied for many years for 
example in Refs. [21-34] by employing the isobar model and/or other rather different theoretical 
approaches, but it is still in the early stage for both the theoretical studies and the experimental 
measurements of such kinds of three-body decays. For example, the factorization for such three- 
body decays has not been verified yet, and many important issues remain to be resolved. 

In Ref. [25], for instance, the authors studied the decays of 5 —?■ Kinr by assuming the va¬ 
lidity of factorization for the quasi-two-body B —)■ {KTr)s,pTi Kinr and introducing the scalar 
and vector form factors to describe the matrix element < J’f“7r+|(sc/)y_^|0 >. 

From the viewpoint of the authors of Ref. [26], a suitable scalar form factor could be developed 
by the chiral dynamics of low-energy hadron-hadron interactions, which is rather different from 
the Breit-Wigner form adopted to study B ^ an. 

In Refs. [27-29], the authors calculated the branching ratios and direct CP violation for the 
charmless three-body hadronic decays B 3h with h = (tt, K) by using a simple model based 
on the factorization approach. They evaluated the non-resonant contributions to the considered 
decays in the framework of heavy meson chiral perturbation theory (HMChPT) with some mod¬ 
ifications, while describing the resonant contributions by using the isobar model in terms of the 
usual Breit-Wigner formalism. The strong phase 0, the parameter and the exponential factor 
g-«NRPs (pi-i-Pj) are introduced in their works [29] in order to accommodate the data. 

In PQCD factorization approach, however, we study the three-body hadronic decays of B me¬ 
son by introducing the crucial non-perturbative input of the two-hadron distribution amplitude 
(DA) [35] and use the time-like form factors to parameterize these two-hadron DAs. In our 
opinion, a direct evaluation of hard 6 -quark decay kernels, which contain two virtual gluons at 
leading order (FO), is power-suppressed and not important. When there is at least one pair of light 
mesons having an invariant mass below O(AmR) [21] ( here A = — mb being the B meson 

and h quark mass difference), the contribution from this region is dominant. The configuration 
involves two energetic mesons almost collimating to each other, in which three-body interactions 
are expected to be suppressed. However, the relative importance of the contributions from the 
two hard gluon exchanges and from the configuration with two collimating mesons still depend 
on specific decay channels and kinematic regions considered. It seems reasonable that the dynam¬ 
ics associated with the pair of mesons can be factorized into a two-meson distribution amplitude 
^hih 2 [35]. One can describe the typical PQCD factorization formula for a i? —)■ 6 - 16 - 26^3 decay 
amplitude as the form of [ 21 , 22 ] 

A = ® H ® ^h-i_h2 ® ( 1 ) 

where the hard kernel H describes the dynamics of the strong and electroweak interactions in 
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FIG. 1. Typical Feynman diagrams for the three-body decays —>• r^cvr+vr with q = {d,s), and the 

symbol • denotes the weak vertex. 


three-body hadronic decays in a similar way as the one for the two-body hadronic B —)■ hih 2 
decays, the functions ^hih 2 are the wave functions for the B meson and the final state 

mesons, which absorbs the non-perturbative dynamics in the process. Specifically, ^hih 2 is the 
two-hadron (hi and /t 2 ) DAs proposed for example in Refs. [35-37], which describes the structure 
of the final state hi — h 2 pair, as illustrated explicitly in Fig. 1. 

By employing the PQCD approach, the authors of Ref. [38] studied the —)■ 7r^7r”'"7r" and 
A^7r+7r“ decays, evaluated the direct CP asymmetries by fitting the time-like form factors and 
the rescattering phases contained in the two-pion distribution amplitudes to relevant experimental 
data, the resulted PQCD predictions agree well with the LHCb measurements [2, 3]. In Ref. [38], 
however, only the non-resonant contributions to the time-like form factors were taken into account, 
the regions involving intermediate resonances are not considered. In the new work [39], by 
parameterizing the complex time-like form factors which include both resonant and non-resonant 
contributions, the authors studied the three-body decays Bg —)■ J/'0/o(98O) [/o(980) — )■ tt+tt”] 
and Bg —)■ /o(980) [/o(980) —)■ 7r+7r“]/i+yU“ decays by using the S'-wave two-pion DAs. 

In recent years, significant improvements for understanding the heavy quarkonium production 
mechanism have been achieved [40]. The meson rjc and J/ifj have same quark content but with 
different spin angular momentum. Following Ref. [39], we here will study the three-body hadronic 
decays 73°^^ —)■ rjcfo —)■ r]c[fo tt+tt”]. We will consider the S'-wave resonant contributions to 

the decay B^ —)■ ?7c/o(500) — )■ ?7c(7r+7r“), as well as the decays Bg — )■ rjcfoiX) — )■ ?7c(7r+7r“) with 
/o(X) = (/o(980),/o(1500),/o(1790)). Apart from the leading-order factorizable contributions, 
we also take into account the NLO vertex corrections to the Wilson coefficients. In Sec. II, we 
give a brief introduction for the theoretical framework and present the expressions of the decay 
amplitudes. The numerical values, some discussions and the conclusions will be given in last two 
sections. 


II. THE THEORETICAL FRAMEWORK 

By introducing the two-pion DAs, the B^^^ —)■ ricTi^'K~ decays can proceed mainly via quasi- 
two-body channels which contain scalar or vector resonant states as argued in Refs. [21, 38]. We 
firstly derive the PQCD factorization formulas for the 73°^^ —>■ ric'K^'K~ decays with the inputs of 
the S'-wave two-pion DAs. We made an hypothesis that the leading-order hard kernel for three- 
body 73 meson decays contain only one hard gluon exchange as depicted in Fig. 1, where the 73° 
or 73° meson transits into a pair of the vr^ and tt^ mesons through an intermediate resonance. The 
Figs. 1(a) and 1(b) represent the factorizable contributions, while the Figs. 1(c) and 1(d) denote 
the spectator contributions. 

In the light-cone coordinates, we assume that the light final state “two pions” and rjc is moving 
along the direction of = (1, 0, Ot) and n_ = (0,1, Ot), respectively. The 73°^^ meson momen¬ 
tum pb, the total momentum of the two pions, p = Pi+ P2, and the pc momentum p^ are chosen 
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as 


Pb — Ot), 


ms 


^r;,0T), p3 = ^(r^r;,0T), 


( 2 ) 


where denotes the meson mass, the variable r] is defined as rj = a;^/[(l — r^)m^] with 
the mass ratio r = rUri^/mB, the variable fj = 1 — rj and the invariant mass squared = 

m^( 7 r+ 7 r“) of the pion pair. As shown in Fig. 1(a), the momentum of the spectator quark in the 
B meson, the momentum k = zp^ and ^3 = are of the form of 


— ( 0, -^Xb, kBT 


k = zp'^ = 


rriB 

V2 


z{l -r 2 ), 0 ,/cT , 


h = X 3 P 3 = 


tub 2 B 2 b - , 

-^r X 3 , ^3T 


where the momentum fraction xb, z and X 3 run between zero and unity. 
The wave function of B/Bg meson can be written as [14, 15, 41] 


( 3 ) 


= 


+ mB)750B(ki). 


( 4 ) 


Here we adopt the B-meson distribution amplitude 0 b( x, b) in the PQCD approach widely used 
since 2001 [14, 15,41] 


(l>B{x,b) = Nbx‘^{1 


a;)^exp 


■ 

2a;| 


^i^B bf 


( 5 ) 


where the normalization factor Nb depends on the value of cjb and /b and defined through the 
normalization relation dx (pBix,b = 0 ) = /b/(2 \/6). ub is a free parameter and we take 
ub = 0.40 ± 0.04 GeV and ub^ = 0.50 ± 0.05 GeV in the numerical calculations. 

For the pseudoscalar meson pc, its wave function can be written as 


4^ 






( 6 ) 


here the twist-2 distribution amplitude 0^ and the twist-3 distribution amplitude pg take the form 
of [42] 


0^(a:) 

i>g{x) 


9.58 

1.97 




2y/W, 

U 


:a;(l — x) 


x{l — x) 


H 0.7 


_l- 2 . 8 x(l -x) 


a:(l — x) 


-[ 0.7 


2y/Wl [1 - 2 . 8 a:(l -x) 


where is the decay constant of Pc meson. 

The S'-wave two-pion distribution amplitude have been defined in Ref. [43] 


( 7 ) 


<|)^ = 


v/2iV; 


7=0 

vu 
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where ( = pf /p~^ is the momentum fraction of the in the pion pair, the asymptotic forms of 
the individual DAs in Eq. (8) have been parameterized as [35-37] 


$ 


7=0 

vu=— 


s 

tv=+ 


9FJ^ 

Fs{uj^) 


02 ^Z{\ 


Fsjuj^) 

2y/M'c 


(l-2z), 


z)(l-2z), 


( 9 ) 


with the time-like scalar form factor Fsiu:^) and the Gegenbauer coefficient 02 “°. For simplicity, 
we here denote the distribution amplitudes ce^), C, and C, by 

00 , 0s and 0t, respectively. 

Following the FHCb collaboration [5, 6, 8] \ we also introduce the S-wave resonances into the 
parametrization of the function Fs{u'^), so that both resonant and non-resonant contributions are 
included into the S'-wave two-pion wave function . For the ss component in the Bg —)■ 
decay, we take into account the contributions from the intermediate resonant /o(980), /o(1500) and 
/o(1790) as in Ref. [8]. We use the Flatte model [44] for /o(980) as given in Eq. [18] of Ref. [8], 
and the Breit-Wigner model for /o(1500) and /o(1790). The ss component of the time-like scalar 
form factor , consequently, can be written as the form of 


^/o(980) ^ ^'^/o(980) (ffTTTrPTrTr T QKKPKK^ 

^ (1500)6'^^ _ 

^/o(1500) ~ ~ *^/o(1500)r/o(1500) 

^_ g3rn^„(i790)e^"^ _ 

"^/o(1790) - '*"i/o(i79o)rjo(i79o)(a;2) ’ 


( 10 ) 


here the three terms describe the contributions from /o(980), /o(1500), and /o(1790), respec¬ 
tively. All relevant parameters in above equation are the same as those being defined previously in 
Refs. [8, 39, 45], such as 


"i/o( 980 ) = 0.97GeV, = 0.167, qkk = 

^/o(i 50 o) = l.SOGeV, mj(,(i 79 o) = 1.81GeV. (11) 

We assume that the energy-dependent width r5(a;^) for a S'-wave resonance decaying into two 
pions is parameterized in the same way as in Ref. [46] 



ms ^ - 4:ml \ ^ ^2 


( 12 ) 


with the pion mass = 0.13 GeV, the constant width Fg with Fs = 0.12, 0.32 GeV for /o(1500) 
and /o(1790) respectively, and the Blatt-Weisskopf barrier factor Fr = 1 in this case [6]. 

For the dd component Ff(a;2), only the resonance /o(500) or the so-called a meson in liter¬ 
ature is relevant. Because the resonance /o(500) is complicated and has a wide width, we here 


* In their analysis [5, 6, 8], the LHCb collaboration used the Flatte model [44] for the description of /o(980), the 
Breit-Wigner model for /o(1500) and /o(1790). 
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parameterize the /o(500) eontribution to the sealar form faetor for the dd eomponent in two differ¬ 
ent ways: the Breit-Wigner and the D. V. Bugg model [47], respeetively. Following Refs. [5, 45], 
we firstly adopt the Breit-Wigner model with the pole mass mj(,( 5 oo) = 0.50 GeV and the width 
r/o( 5 oo) = 0.40 GeV, 



_^^/o(500)_ 

’^/o(500) ~ ~ *^/o(500)r/Q(500) 


(13) 


The parameters c, c* and 0* with i = (1, 2, 3) appeared in Eqs. (10) and (13) have been extraeted 
from the LHCb data [ 8 ], 


Cl = 0.900, C 2 = 0.106, C 3 = 0.066, c = 3.500, 

01 = -^, 02 = ^, 03 = 0. (14) 

Seeondly, we parameterize the form faetor of /o(500) with the D. V. Bugg resonant lineshape [47] 
in the same way as in Ref. [48] 


Tn(s) = Mri(s) 




2 S ~ SA 


M2 - 




2=1 


-1 


(15) 


where s = = m^{7i+7i ), ji(s) = i 2 + pi In 

relevant funetions in Eq. (15 ) are the following 


, the funetions gf{s),Ti{s) and other 


gl{s) = M{hi + h 2 s) exp[-(s - M^)/A], 

M ri(s) = ^i(g) ^2 

M 12 ( 5 ) = Q.<6g‘l{s){s/M"^) exp(-a|s - 4m|^|)p2(s), 

M 13 ( 5 ) = Q.2gl{s){s/M'^) exp(-a|s - 4mJ|)p3(s), 

M r 4 (s) = M g4T, p4n{s)/ p 4 n{M‘^), with p4Tr{s) = 1.0/[1 + exp(7.082 — 2.845s)], (16) 


For the parameters in Eqs. (15,16), we use their values as given in the fourth eolumn of Table I in 
Ref. [47]: 

M = 0.953GeV, s^ = 0.41 m^, hi = 1.302GeV2, 

b 2 = 0.340, A = 2.426GeV^ ^ 4 ^ = O.OllGeV. (17) 


And the parameters pi, 2,3 in Eq. (16) are the phase-spaee faetors of the deeay ehannels tttt, KK 
and pp respeetively, and have been defined as [47] 





(18) 


with rrii = m^r, m 2 = mx and m 3 = m^. 

The differential deeay rate for the 5°^^ —)■ pcTi^Ti~ deeay ean be written as [45] 


^ ^ ^I:PiIIp3| 14,2 
duj ^4(27r)3mg 


( 19 ) 
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where oj = m(7r+7r ), \pi\ and Ipal denote the magnitudes of the 7r+ and rj^ momenta in the 
eenter-of-mass frame of the pion pair, 


|pi| = -4m2±, 


The deeay amplitude for the deeay rjcTT^'K~ is of the form 

^ r^.TT+TT-) = 

- VtlVtdits) {F'^^ + F^^ + , 


( 20 ) 


( 21 ) 


where the funetions and and ) denote the amplitudes for the 

B/Bg meson transition into two pions as illustrated by Fig. 1(a) and 1(b) ( Fig. 1(e) and 1(d)): 

pi poo 

F^^ = ^T[Cpm%fn^ / dxsdz / bs dbs b db (j)B{xB,bB) 


X 


\/r]{l - r2) [[(1 - 2z)rj + r^(l + 2zr])](j)s - (r^(l - 2zr]) + {2z - l)r])4>t\ 


[ — (1 + z)r] + r^(l + 2zr] — 27]) — r‘^zri\ 0o 


ai(fa) Ee{ta) ha{x B, Z ,b B ,b) 


2^Jr]{l - r2) [r^{r] - xb) + (1 - - (1 - r'^) [w + “ XB)](t>o 


X ai(4) Ee{tb) hb{xB,z,bB,b) L 


jpfLL _ ^LL\ 


l<2.1—)-a25 


= —P 


LL\ 


|ai—)-a3 1 


( 22 ) 

(23) 


X 


= 327 rCi?m ^/\/6 / dxBdzdx^ / bBdbB b^db^ (j)B{xB,bB) 

Jo Jo 

^^(1 _ ^2) ^^2^^ ('l _ ^^('X _ 

- [2r^(l - Xa)?) + (1 - r^)2:?) - - ArrpipgCt)^ 

+ (r) + r^) [(1 - r^)(l - X3 - Xb) + 3:3(1 - 2r^)r7 - (1 - r^)(l - z)p + 

- rrc(l - + ? 7 )V:s 0 o <^2(4) EniF) H^xb, z, X3, bB, ^3) 

^77(1 - r2) [ 2 r^X 3 r 7 - + -^(1 - r'^)fi\'ijjv(l)s + [-2(1 - r‘^)f] + r‘^XB]jJv 4 ’t 


- [r] + r‘^{r]- p)] [x^p - xb + + z{l - r^)]'0«0o 


X C2{td) En{td) hd{xB,Z,X3,bB,b3) [>, 
= M^^\c2^C4+Cio^ 


(24) 


( 25 ) 
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= 32 TTCFm%/V g I dxB dz dx^ I bBdbB b^db^ 0s(a^s,&s) 

Jo Jo 

x| ^77(1 - r2) [ 2 r'^{l - X 3 )f] + z{l - r^)!] - r‘^XB]'ilJv(l)s - 4 :rrc'ilJs(l>s 
+ \/r]{l - r2) [2;(1 - r'^)f] + 

(-r^ + r/ + 2 r‘^r]) [(1 - a;3)r^ + z{l - r^) + (1 - x^)?] - XB]i>v(j)o 


rrc(l + r] — r‘^)'ijjs 


[Coitc) + Csitc)] Er,{tc) K{xb, Z, X3, bB, 63) 


^ 77(1 - r 2 ) [ 2 r^X 3 r 7 + (1 - r‘^)z7] - r‘^XB]i>v(l>s - [- 2(1 - r‘^)r] + r‘^XB\i>v(t>t 


— (r^ + 77 ) (1 — r‘^){x 3 rj — Xb) + X 3 r‘^rj + ? 7 z(l — r^) 7 /;, 
X [C'6(td) +C8{td)] En{td) hd{xB,Z,X3,bB,b3) k 


( 26 ) 


where Cf = 4/3, rc = rric/mB, and a* are the combinations of the Wilson coefficients Cf. 

_ , 6^2 _ /-f I I ^4 + C'lo _ ^ ^ , C'e + C's 

ai — Ci + —, 02 — W 3 + wg H-^-, 03 — 05 + (77 H ^-. (27) 

The explicit expressions of the evolution factors (-Eg(4), Ee{h), En{tc), En{td)), the hard func¬ 

tions {ha, hb, he, hd) and the hard scales {ta, h, tc, td), appeared in Eqs. (22-26), can be found in 
the appendix of Ref. [39]. 

For the factorizable emission diagrams Fig. 1(a) and 1(b), the NFO vertex corrections can be 
taken into account through the inclusion of additional terms to the Wilson coefficients [11, 12, 49, 
50]. After the inclusion of the NFO vertex corrections, the Wilson coefficients Oi, 02 and 03 as 
defined in Fq. (27) will be modified into the following form 




0'2{h) — 6^3(/^) + Cg{ix) -f 

2 as(/i) 


Svr 

C^{li) + (7io(/r) 


61 n—— 9 - 1 --^ [ dx'ijjy{x)g{x) 

h trie Jo 


X U + 


Stt 


61 n—— 9 - 1 --^ [ dx'ipy{x)g{x) 

h trie Jo 


O'oih) — + 6^7 (/^) + 

2as{g) 


Ceih) + 


X 1 


Svr 


61 n—— 3 - 1 -^^ [ dx'ipy{x)g{l — x) 

h tve Jo 


( 28 ) 


( 29 ) 


( 30 ) 


Since the emitted meson 77 ^ is heavy, the terms proportional to the factor 2 ; = can not 

be neglected. One therefore should use the hard-scattering functions g{x) as given in Ref. [51] 
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instead of the one in Ref. [49], 


9{x) = 


3(1 - 2x) 

1 — X 

x"^ {\n.[l — z] — in) 2 


(1 — (1 — x)z)'^ 


ln[x] + 3 [ln(l — z) — in] 
+ xz'^ \n[xz\ 


2z{l — x) 2 a: 2 ;(ln[l — z] — in) 


1 — zx 

X 


1 — (1 — x)z 
1 — X 


(1 — (1 — x)z)'^ (1 — xz)'^ 


2 xz \yi[xz\ 


1 — (1 — x)z 1 — xz 


(31) 


where 2 ; = For Bg —>■ rjcn^n deeay, the mass in above equations should be 

replaeed by the mass ms^. 


III. NUMERICAL RESULTS 

In numerieal ealeulations, besides the quantities speeified before, the following input parame¬ 
ters (the masses, deeay eonstants and QCD seale are in units of GeV ) will be used [38, 45] 

= 0.326, mfio = 5.280, = 5.367, = 2.9836, 

= 0.140, rriT^o = 0.135, = 4.8, = 1.275, = 0.095, 

fs = 0.19 ± 0.02, fs, = 0.236, r^o = 1.519 ps, tb, = 1.512 ps. (32) 

The values of the Wolfenstein parameters are the same as given in Ref. [45]: A = 0.814+0 024 , A = 
0.22537 ±0.00061, p = 0.117 ±0.021, f) = 0.353 ±0.013. For the Gegenbauer eoeffieient we use 

4=^ = 0 . 2 . 

In Fig. 2(a), we show the eontributions to the differential deeay rate dB{Bs —)■ ? 7 cVr+ 7 r“)/da; 
from eaeh resonanee /o(980) (the blue solid eurve), /o(1500)(the red solid eurve) and /o(1790) 
(the dots eurve), as a funetion of the pion-pair invariant mass u = m{n^n~)- For the eonsidered 
Bg deeay, the allowed region of oj is 4m^ < < {Mb^ — In Fig. 2(b), furthermore, we 

show the eontribution to the differential deeay rate dB{B^ n~)/du from the resonanee 

/o(500) as a funetion of m{n^n~) too, where the red and blue line shows the predietion obtained 
by using the Breit-Wigner model and the D. V. Bugg model respeetively. For B —)■ ? 7 c 7 r+ 7 r” deeay, 
the dynamieal limit on the value of oj is 4m^ < < {Mb — 

For the deeay B —Pc/o(500) — ric'^^n~, the PQCD predietion for its branehing ratio with 
Breit-Wigner form is 

B(B<‘ ^ ,,/„(500)[,r+,r-]) = X Ifl-'. (33) 

When we use the method of D. V. Bugg, the PQCD predietion for its branehing ratio is of the form 

B{B^ ^ pJo(500)[7r+7r-]) = [2.3i;°;i(a;s)l[i:S[i(ar°)l°;°?(m,)] x 10-^ (34) 

where the three major errors are indueed by the uneertainties of ojb = (0.40 ± 0.04) GeV, 4^^ = 
0.2 ± 0.2 and rric = (1.275 ± 0.025) GeV, respeetively. 

For the deeay mode Bg —)■ pcfoi^) Vc{t^^'^~)s, when the eontribution from eaeh resonanee 
/o(980), /o(1500) and /o(1790) are ineluded respeetively, the PQCD predietions for the branehing 
ratios for eaeh ease are the following, 

B{Bg Pc/o(980)[7r+7r“]) = [3.37l|(;?^(a;Bj;[(;|^^(a2=°)lo;oo(mc)] x 10■^ (35) 

B{Bg ^ pJo(1500)[7r+7r-]) = [6.76l};^?(a;Bj;no(«2=°);o:ooK)] x lO"®, (36) 

B{Bg ^ pJo(1790)[7r+7r-]) = [l.97l°i9(u;Bj;E5:^J(a^=°)l[J:|J}(me)] x 10"®, (37) 
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FIG. 2. The m(7r+7r“)-dependence of the differential decay rates dB/dco for (a) the contribution from the 
resonance /o(980), /o(1500) and /o(1790) to —)• rycvr+vr" decay, and (b) the contribution from /o(500) 

to —)• decay in the BW model ( red curve) or the Bugg model ( blue curve) . 


where the three major errors are indueed by the uncertainties of = (0.50 ± 0.05) GeV, = 
0.2 ± 0.2 and nic = (1.275 ± 0.025) GeV, respectively. The errors induced by the variations of the 
Wolfenstein parameters and other inputs are very small and have been neglected. If we take into 
account the interference between different scalars fo{X), we found the total branching ratio: 

B(B. ^ qMX)l,T+7r-]} = [5.02i;;ig(ws.)«;“(of x lO'^. (38) 

The interference between /o(980) and /o(1500), /o(980) and /o(1790), as well as /o(1500) and 
/o(1790), will provide a contribution of 3.29 x 10“®, 5.59 x 10“® and —1.16 x 10“® to the total 
decay rate, respectively. 

From the curves in Fig. 2 and the PQCD predictions for the decay rates as given in Eqs. (34-38), 
one can see the following points: 

(i) For Bg 7]cfo{X) ric{fo{X) tt+tt") decay, as illustrated clearly by Fig. 2(a), the 
contribution from the resonance /o(980) is dominant (~ 70%), while the contribution from 
/o(1790) is very small (~ 4% only). The interference between /o(980) and /o(1500), as 
well as /o(980) and /o(1790), are constructive and can provide ~ 20% enhancement to the 
total decay rate. The interference between /o(1500) and /o(1790), however, is destructive, 
but very small (less than —2%) in size. 

(ii) For B — > ?)c/o(500) — )■ r^cvr+vr” decay, the PQCD predictions for its branching ratios are 
around 2 x 10“® in magnitude when we use the Breit-Wigner or the D. V. Bugg model 
to parameterize the wide /o(500) meson. The model-dependence of the differential decay 
rate dB/duj{B —)■ r7c/o(500)(7r+7r“)), as illustrated in Fig. 2(b), are indeed not significant. 
Although the central value of PQCD predictions based on the D. V. Bugg model are mod¬ 
erately larger than the one from the Breit-Wigner model, but they are still consistent within 
errors. 

(iii) The B/Bg ^ r]cfo{X) rjc(7r^7r~) decays are similar in nature with the decays B/Bg ^ 

J/ihfniX) J/'0(vr+7r“) studied previously in Ref. [391. We find numerically BiB 
r^cvr+vr-) : B{B ^ 0.3 : 1. 
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IV. SUMMARY 

In this paper, we studied the eontributions from the S'-wave resonant states fo{X) to the 
r]c7i^7T~ deeays by employing the PQCD faetorization approaeh. We ealeulated the differential 
deeay rates and the branehing ratios of the deeay —)■ rjcfoi^OO) —)■ ric{n~^Tr~), the deeays 
Bs Vcfo{X) r]^{7T+7T-) with fo{X) = /o(980),/o(1500) and /o(1790) respeetively. By 
using the S'-wave two-pion wave funetion the resonant and non-resonant eontributions to the 
eonsidered deeays are taken into aeeount. The NLO vertex eorreetions are also ineluded through 
the redefinition of the relevant Wilson eoeffieients. 

From analytieal analysis and numerieal ealeulations we found the following points: 

(i) For the branehing ratios, we found 

^ rjJoibOO) ^ r^eTT+TT-) = x 10-^ (39) 

l^Bugg^^O ^ ^J^(500) ^ r^^VT+TT-) = (2.311^8) ^ 10-^ (40) 

BiB, ^ VcfoiX) ^ r^.TT+TT-]) = (5.02l|;«) x 10-^ (41) 

where the individual errors have been added in quadrature. For the deeay rate B{Bs 
ric'K^'K~), the eontribution from the resonanee /o(980) is dominant. 

(ii) For B —)■ rjcfo{500) decay, we used the Breit-Wigner and the D. V. Bugg model 

to parameterize the wide /o(500) meson respectively but found that the model-dependence 
of the PQCD predictions are not significant. 

(iii) The considered decays with the branching ratio at the order of 10“® ~ 10“^ could be mea¬ 
sured at the ongoing LHCb experiment. The formalism of two-hadron distribution ampli¬ 
tudes, consequently, could be tested by such experiments. 
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